CU-TP-718 



Excited states in the full QCD hadron spectrum on a 16^ x 40 lattice 

Dong Chen'' * 

''Department of Physics, Columbia University, New York, NY 10027, USA 



in 

o 
O 

m 

m 
in 
o 
o 

in 



We report the hadron mass spectrum obtained on a 16^ x 40 lattice at /3 = 5.7 using two flavors of staggered 
fermions with ma = 0.01. We calculate the masses of excited states that have the same quantum numbers as the 
TT, p and N. They are obtained by a combined analysis of the hadron correlators from sources of size 16^ and 8^ 
We also report on the hadron spectrum for a wide range of valence quark masses. 



1. INTRODUCTION 

The hadron mass spectrum is of fundamen- 
tal interest in lattice QCD. The Columbia group 
has done hadron mass calculations with two fla- 
vors of staggered dynamical quarks for ma = 
0.01,0.015,0.02,0.025 on a 16^ x 32 volume and 
ma = 0.01,0.025 on 32^ both at ,5 = 5.7 [1- 
3]. We have also investigated ma = 0.004 on a 
16^ X 32 volume at /3 = 5.48 [4]. 

At Lattice '94, we reported a new studj^ on a 
16^ X 40 lattice with ma = 0.01 at /3 = 5.7 [5]. 
Since then, we have extended our calculation to 
4900 micro-canonical time units. We systemati- 
cally study the effects of source size on the deter- 
mination of the mass values. We obtain the ex- 
cited states by combined analysis of hadron prop- 
agators from different sized sources. We also re- 
port on the hadron mass spectrum with a variety 
of valence quark masses and compare our results 
with our earlier runs. 

2. SIMULATION 

Table 1 hsts the parameters of our calculation. 
We used the 'R'-algorithm of Gottlieb, et. al. [6] 
for our evolution, along with their notation. W^e 
have collected 4900 micro-canonical time units 
during a 7.5 month run on the 256-node Columbia 
parallel computer. Hadron masses were measured 
every 6 time units. For a single gauge configura- 



tion, we did 5 measurements with the source set 
on 5 different time slices and averaged the result 
over time shces (AOTS). 

Table 1 

Simulation parameters 



volume 

a 

dynamical^ 

length 

trajectory length 
step size 

CG stopping condition 
avg CG steps 
minutes / trajectory 

hadron measured every 
hadron source types 

hadron source sizes 
valence quark masses 



16=* X 40 
5.7 
0.01 
250-1-4650 

0.5 
0.0078125 
1.13 X 10-^ 
660 
12.5 

6 

wall, 2Z 
wall, Z 

l6^l2^8^4^ point 

0.004, 0.01, 0.015, 0.02 
0.025, 0.05, 0.07 



*This work was done in collaboration with Shailesh Chan- 
drasekharan, Norman H. Christ, Weonjong Lee, Robert D. 
Mawhinney and Decai Zhu and was supported in part by 
the Department of Energy. 



In our hadron mass calculations, the 2Z source 
is set on a fixed t slice and is non-zero when only 
all {x,y,z) are even; our Z source is non-zero for 
all {x,y,z). Wc set our 2Z and Z wall sources 
across the entire spatial volume (16^) as well as 
on smaller cubic spatial volumes (12^, 8^, . . .). We 
also calculate quark propagators using different 
valence quark masses, as hsted in Table 1 . Out 
of all the possible combinations of source types, 
sizes and valence quark masses, we have a total 
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of 19 combinations for mesons, 21 combinations 
for the nucleon and 10 combinations for the A. 

Figure 1 shows the evolution of (xx) during 
our run. Discarding the first 250 time units for 
thermahzation gives (xx) = 0.0274(2). This is 
consistent with the vakie 0.0277(3) from our ear- 
Uer 16^ X 32 run [1,2]. 
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Figure 1. The evolution of the chiral condensate. 



3. RESULTS 

We reported source size effects on the effective 
mass plots for hadrons in ref [5]. Here we repro- 
duce two such plots with more statistics. Fig- 
ures 2 and 3 show from bottom to top the effec- 
tive mass plots of the tt, p and N for 16'^ and 8^ 
sources. These figures demonstrate that 16^ is the 
optimal source size for this calculation. However, 
they also show that close to the source, there are 
signals for excited states since the effective masses 
are heavy. 

Since the excited states have the same quantum 
numbers as the lowest states, a straight-forward 
fit to both lowest and excited states tends to have 
oscillations between the mass parameters result- 
ing in poorly determined masses. We thus make 
combined fits using both 16^ and 8^ sources. The 
16^ data are only fitted to the lowest states while 
the 8^ data are fitted to both lowest and excited 
states with the lowest states mass set identical to 
those for the 16^ data. (Each particle comes with 
its opposite parity state for staggered fermions. 
The tt's parity partner does not exist in the nu- 
merical data.) As an example, for the pion we 
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Figure 2. Tr,p,N effective masses for 16^ source. 
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Figure 3. Tr,p,N effective masses for 8^ source. 

perform a one state fit to the 16^ source with 
P^it) = Aocosh{moit-Nt/2)), (1) 
and a two state fit to the 8^ source with 



P^{t) = A'oCosh{mo{t-Nt/2)) 
+ A'^ coshim'^ (t - Nt/ 2)). 



(2) 



The two mo's are set to be the same; thus the 
combined fit has 2-1-4 — 1 = 5 parameters. For 
all other particles, we make a 2 state fit to the 
16^ source and a 4 state fit to the 8^ source with 
the two lowest mass values identical and hence 
the combined fit has 4 + 8 — 2 = 10 parameters. 
The degree of freedom (dof) for our fits is thus 
2 X {Nt/2 — tmin) minus the number of fitting pa- 
rameters. We include the full correlation matrix 
in ah our fits. 

Figures 4 and 5 show the combined fitting re- 
sults for the TT and p. These are not effective mass 
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Figure 4. The tt and its excited state from a com- 
bined fit for the 16^ and 8^ sources. 
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Figure 6. The Ai effective masses in Coulomb 
gauge (o) and Landau gauge (x) for a 16^ source. 
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Figure 5. The p and its excited state from a com- 
bined fit for the 16^ and 8^ sources. 



plots, but the fitted masses from a distance tmin 
up to Nt/2 as a function of tmin- There are clear 
signals for excited states at small tmin- 

Table 2 is a Hst of hadron masses for our cur- 
rent 16^ X 40 calculation, as a function of the 
valence quark masses we used to calculate our 
hadron propagators. All our fits are correlated 
fits from tmin to Nt/2 with errors determined 
from the jackknife method. The pion mass is ob- 
tained from a single state fit with the 2Z source. 
All other mesons and the nucleon come from two 
state fits with the 2Z source. Our A's come from 
two state fits with the Z source. Ai and Aq are 
two equivalent lattice operators corresponding to 
the two A2 operators in ref [7], respectively. All 
our sinks are local except for the A's which are 



necessarily non-local in the staggered fermion for- 
malism. All our source and sink time slices are 
fixed to Coulomb gauge. For the A, we also mea- 
sure in Landau gauge for comparison. 

Figure 6 is an effective mass plot for the Ai. 
The Aq plot is very similar. For small t, the A 
masses measured in Landau gauge are quite dif- 
ferent from those measured in Coulomb gauge. 
They are about 3 standard deviations lower and 
are in closer agreement with the nucleon. This 
raises a question about the validity of using Lan- 
dau gauge, non-local operators in general, since 
the gauge-dependent correlation between differ- 
ent time slices in Landau gauge may alter the 
signal one wants to obtain. 

Table 3 lists the excited states in our calcu- 
lation. We only measured 8^ source data for 
ma = 0.01. The tmin used is less than in table 2 
so that we have a better estimate for the excited 
states. Were we to use the same tmin in table 2 
as in table 3, we would have < 1% change in the 
mass values but a larger x^- The lowest states 
from our combined 16'' source and 8'^ source fits 
also agree with those in table 2 to 1%. All our fits 
have x^/dof < 2. There is no apparent excited 
state for the 71-2 in our data. 

Compared with our earlier 16^ x 32 calcula- 
tion [2], where we have in„a = 0.252(3), m^a = 
0.454(4) and rriNa = 0.692(6), our current 16'^ x 
40 mass spectrum is lighter. The nucleon is es- 
pecially light, about 3 standard deviations lower 
than the earlier result. In our current calculation, 
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Table 2 

Hadron mass results: 16^ x 40,/3 = 5.7, ma = 0.01 







0.004 


0.01 


0.015 


0.02 


0.025 


0.05 




0.07 






















m„a 


11 


0.167(2) 


0.248(2) 


0.299(2) 


0.343(1) 


0.382(1) 




15 


0.643(1) 




8 


0.210(7) 


0.285(3) 


0.335(3) 


0.380(3) 


0.421(3) 




11 


0.702(2) 






0.320(13) 


0.424(11) 


0.474(10) 


0.516(9) 


0.553(9) 






0.811(4) 


mpCi 


6 


0.407(5) 


0.438(4) 


0.465(4) 


0.492(4) 


0.519(4) 




12 


0.754(3) 






0.590(36) 


0.599(11) 


0.626(9) 


0.654(9) 


0.683(9) 






0.940(29) 


nip^a 


6 


0.398(10) 


0.433(6) 


0.463(5) 


0.492(4) 


0.521(4) 




10 


0.758(3) 


rua^a 




0.513(13) 


0.564(11) 


0.598(11) 


0.629(10) 


0.660(10) 






0.905(9) 


niNa 


6 


0.600(10) 


0.661(7) 


0.708(6) 


0.753(6) 


0.797(5) 




11 


1.163(5) 


niNia 




0.674(18) 


0.757(12) 


0.815(11) 


0.866(12) 


0.914(11) 






1.256(13) 




6 


0.734(33) 


0.746(16) 






0.852(8) 


1.040(6) 








6 


0.752(47) 


0.729(14) 






0.843(9) 


1.035(6) 







Table 3 

Excited states. 









X' 




0.68(4) 


7 


49(17) 












0.73(6) 


6 


18(10) 


nipa 


0.90(3) 


4 


36(19) 




0.98(12) 








0.91(7) 


4 


40(16) 




0.94(5) 






niNCi 


1.17(2) 


4 


32(17) 


rriNia 


1.19(3) 







The excited states are labeled as their correspond- 
ing lowest states. Other particles that have differ- 
ent quantum numbers in the continuum can also 
have projections onto these states on the lattice. 



the masses do not depend on the upper bound 
of the fitting range within errors. The length in 
micro-canonical time units where we do hadron 
measurements for the new run is about 3 times 
longer than the old one, although the old run has 
more averages for a single gauge configuration. 
We would argue that in our previous run, the er- 
rors are underestimated because our simulation 
time is not sufficiently large to accurately com- 
pute the effects of autocorrelations. 



Using the p mass of 770 MeV to fix the scale, 
we find at ma = 0.01 that m^ = 436(4) MeV, 
mjv = 1162(12) MeV and m^ = 1297(28) MeV 
(averaging over our Coulomb gauge Aq and Ai). 
The excited states have = 1195(70) MeV, 
= 1582(53) MeV and mjv = 2057(35) MeV. Our 
lattice spacing is a~^ = 1757(16) MeV and the 
lattice size is La = 1.80(2) fm. 
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